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ABSTRACT 



We have analysed temporal evolution in centroid fre- 
quencies and splitting coefficients of solar /-modes 
obtained from MDI/SOHO. The data were divided 
into 20 sets covering a period from May 1, 1996 to 
August 31, 2000. The variation in frequencies is es- 
timated to be 68 nHz over the period of four years 
which includes the rapidly rising phase of the solar 
cycle 23. This change is much smaller than that ob- 
served for p-mode frequencies. It is also noticed that 
the /-mode frequencies appear to be weakly corre- 
lated with solar activity indices as compared to the 
p-mode frequencies. We have also inferred the rel- 
ative change in the solar radii and notice a 1 year 
periodicity which may be associated with the solar 
cycle variation. 
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1. INTRODUCTION 



With the availability of precise frequencies from both 
Global Oscillation Network Group (GONG) and So- 
lar Heliospheric Observatory (SOHO) instruments, 
refined analysis of solar cycle dependent parameters 
are being carried out. In these studies, p-mode fre- 
quencies have been widely used since these modes 
penetrate deep inside the sun and thus provide a way 
to infer changes beneath the surface. In contrast, 
/-modes, which are the surface gravity modes, are 
less well studied. Most of the efforts in /-mode stud- 
ies have been concentrated in estimating the seismic 
radius of the sun since these mode frequencies are 
indep endent of the stratification in the solar inte- 

pOOOl ; 



rior (pziembowski et al., 1998 



Antia et al. 



OOOaj). Frequency splittings obtained from /-modes 



rave also been used to study the zonal fl ows which 



arc associated with torsional oscillations (Howard 



LaBontc, 198C ) as seen in BBSO data (Woodard & 
bbrecht, 1993). These sub-surface flows were first 



noticed by Kosovichev fc Schou (1997 ) in the Michel- 
son Doppler Imager (MDI) /-m ode splittings . This 



zonal flows towards the equator. Toomre et al. (2000| ) 



analysis was later extended by Schou (1999 ) for a 
longer time period which showed average drift of the 



have confirmed these zonal flows using both p- and 
/-mode splittings derived from MDI data. 

Clearly, we are still lacking a basic understanding 
as to what causes the frequency changes. Keeping 
this in view, we analyse the /-mode frequencies and 
splitting coefficients obtained from the MDI medium- 
£ program for the period which includes the rising 
phase of the solar cycle 23. We look for possible solar 
cycle changes in the mean frequencies and even order 
splitting coefficients. These are further correlated 
with activity indices. We have also estimated the 
variations in the seismic radius for the period covered 
in this study. 



THE DATA 



The results presented here consist of twenty 72 day 
data sets obtained with the MDI instrument on 
board SOHO and covers a period from May 1, 1996 
to August 31, 2000 with two breakdowns in between 
(see Schou, 1999 for more details). Each data set 
consists of centroid frequencies Vn.i and 36 splitting 
coefficients Ofc. For the present work, only /-modes 
are used which samples a region very close to the 
surface. In particular, the analysis is restricted to 41 
common modes in the degree range of 217 < < 286 
and frequency range of 1484 /iHz < ly < 1702 fiHz 
unless otherwise mentioned. 



3. ANALYSIS AND RESULTS 



3.1. Variation of the Sun's Eigen Frequencies 

The variation of i averaged frequency shift as a func- 
tion of frequency for low (May 1 - July 7, 1996) and 
high (June 21 - August 31, 2000) activity period with 
reference to set 6 (April 26 - July 6, 1997) is shown 
in Figure 0. During the minimum activity period, 
dv seems to be independent of v while during high 
activity period, the frequency shift increases sharply 
with frequencies. The mean shift Sv for a given i is 
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Figure 1. The binned frequency shifts for two inde- 
pendent periods: set 1 (May 1 - July 7, 1996) repre- 
sented by diamonds and set 20 (June 21 ~ August 31, 
2000) by stars. The shift has been calculated with ref- 
erence to set 6 (April 26 - July 6, 1997). The error 
bars represent mean error in shift. 



Figure 2. The temporal evolution of frequency shifts 
for p- ( line joined by triangles) and f-modes (line 
connected by diamonds). The shift has been calcu- 
lated with respect to mean of all 20 data sets sepa- 
rately. The dashed line represents the solar activity 
levels denoted here by the International sunspot num- 
ber. For clarity, the error bars are not shown. 



calculated from the relation 



5v{t) 



(t) 
^2 



(1) 



where ai is the error in the frequency measurement. 
The variation of 5v with time is shown in Figure ^. 
Here we notice an oscillatory behaviour with a pe- 
riodicity of about one year. A similar result is also 
found by Antia et al. (2000b| ). The temporal evolu- 
tion of p-mode frequencies (1500 /iHz < v < 3500 
/^Hz) and the activity index represented by the In- 
ternational sunspot number Rj is also shown in the 
same figure. It is evident that the change in p-mode 
frequencies is larger than that of /-mode frequencies. 
It is further clear that the p-mode frequency shifts 
follow the activity index quite closely while the mean 
frequency shifts of /-modes appear to be weakly cor- 
related (also see Table 1). 

To study the correlation between mean frequency 
shifts and solar activity cycle, we have considered the 
foUipwliig dllfoi'iJiiL acLlviLy liidlcus: i?/, Llie IiiLeniu- 
tioifeti- t number obtained from the Solar Gco- 



physical Data (SGD); KPMI, Kitt Peak Magnetic 
Index which represents the line of sight magnetic 
field setrength and averaged over the full disk from 
Kitt Peak magnetograms; MPSI, Magnetic Plage 
Strength Index from Mount Wilson magnetograms 
which represents the absolute values of the magnetic 
field strengths between 10 and 100 gauss; MWSI, 
Mount Wilson Sunspot Index which represents the 
magnetic field strengths above 100 gauss; He I, equiv- 
alent width of He I 10830 A line, averaged over the 
whole disk from Kitt peak and i^io, integrated radio 
flux at 10.7 cm from SGD. Both Pearson's (r^) and 
Spearman's (r^) correlation coefficients are given in 
Table 1. 



Table 1. Correlation statistics for f- and p-mode fre- 
quency shifts. 



Activity 


/-mode 


p-mode 


Index 


rp 


rs 


^ s 


Ri 


0.84 


0.80 


0.98 0.93 


KPMl 


0.88 


0.78 


0.99 0.94 


MPSI 


0.87 


0.77 


0.99 0.95 


MWSI 


0.87 


0.82 


0.91 0.83 


Fio 


0.87 


0.78 


0.99 0.96 


He I 


0.85 


0.79 


0.99 0.97 



It is observed that in case of /-modes, the activity in- 
dices have weak correlation with frequency shifts. In 
comparison, p-modes display higher correlation coef- 
ficients showing stronger correlati on with activity in- 
dices as was earlier poi nted out by Jain, Tripathy, & 
Bhatnagar (2000|) and |Howe, Komm, fc Hill (19991} 



However, it is interesting to note that MWSI with 
other magnetic indices has strong correlation with 
/-modes while it has the weakest correlation with p- 
modes. Thus, it appears that MWSI which repre- 
sents the average magnetic field strenth above 100 
gauss behaves differently with p- and /-modes and 
needs to be investigated more closely as the solar 
cycle progresses. The observation that /-modes in 
comparison with p-modes are weakly correlated with 
magnetic field indices further confirms that magnetic 
fields are present at deeper layers rather than near 
the solar surface. 
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3.2. Seismic Radius 



Since the early nineteen century, the temporal vari- 
ation in the Sun 's diameter has been reported b y 
many observers ( Laclare et al., 1996| ; Noeel, 1997 ). 
The recent measurements report a change of 0.1 arc- 
sec to 1 arcsec in the measured semi-diameter which 
implies a change of 70 to 700 km in the radius. As- 
suming that the change in frequencies is entirely due 
to the change in radius and there is no contribution 
from the magnetic field, convection etc., the average 
frequency difference is related to change in radius by 



1.0 



6u 3 5R 



(2) 



where Sv is the centroid frequencies calculated ear- 
lier and the differences are taken with respect to the 
mean of all 20 data sets. The evolution of the relative 
differences in the seismic radius is shown in Figure 3. 
The maximum change is found to be 29.8 km over 
a period of four years, the seismic radius decreasing 
with increase of activity lavels and thus indicating 
anti-correlation between them. A relative variation 
in the solar radius of the order of 6xl0~^ was first 



reported by Dziembowski et al. (199S) which corre- 
sponds to approximately 4 km over a period of 2 
years around the solar minimum. If this variation is 
correlated with solar activity, we would expect much 
larger change, as reported in this paper, between 
minimum and maxium activity perio d. However, the 
analysis of Dziembowski et al. (2000| ) did not find any 
systematic trend of correlation with suns pot num- 
ber . A si milar study using GONG data (Antia et 
al., |2000aD for the period 1995 to 1998 reported a 
decrease in solar radius by about 5 km, during the 
period in which the solar activity was increasing. It 
has been known that the high wave number modes 
are known to devia .te significantly from the simple 
dispe rsion relation ( Duvall, Kosovichev fc Murawski, 



199^; Antia fc Basu, 1999| ) and the inclusion of higher 
degree in our analysis may cause this large change in 
the solar radius. 

In order to check the stability of our results, we have 
repeated the analysis by considering the frequency 
tables which contains the splitting coefficients only 
up to og and thus lowering the £ values to 119. The 
relative difference in radius in this case is found to 
be 14 km and clearly demonstrates the degree de- 
pendence. We have also estimated the relative dif- 
ference in seismic radius from activity minimum to 
maximum by restricting the common modes between 
100 < £ < 200 and this is shown as dashed line in 
Figure 3. In this case, the change in radius is found 
to be 7 km, consistent with the results of |Antia et 
al 



2000a) which used the same range of £ values 



Since the change in seismic radius reflects the rela- 
tive differences in frequencies, the oscillatory pattern 
visible in Figure 2 persists in Figure 3. However, it 
is not clear if this period is associated with solar cy- 
cle variations since the data used here spans only 4 
years. Also, this period is very close to the orbital 
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Figure 3. Relative differences in the radius of the sun 
as inferred from f-mode frequencies. The solid line is 
the result as derived from frequency tables containing 
splitting coefficients up to a^e and covers the degree 
range of 21 7 to 286. The dashed line shows the result 
using the frequency tables containing splitting coeffi- 
cients up to ae and includes modes of 100 < £ < 200 
only. 
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Figure 4- Average even order splitting coefficients as 
a function of index and time. For each index i, the 
a- coefficients are joined with lines denoting temporal 
evolution. Error bars are 1 a. The averaging was 
done only over common modes. 
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period of the earth and one may expect systematic 
errors of this period being introduced in the data. 



3.3. Even Order Sphtting Coefficients 

The variation of even order coefficients up to aie for 
/-modes are shown in Figure 4. These coefficients 
seem to be uncorrelated with solar cycle except for 
02. In contrast, the even order splitting coefficients 
of p- modes are kno wn to be correlated with solar 



cycle. In particular, Howe, Komm, & Hill (1999) find 



significant correlation between the even a-coefficients 
and the corresponding Legendre coefficients of a fit to 
the Kitt Peak magnetograms. A similar result based 
on the BBSO Ca K plage data for higher orde r was 
recently reported by Dziembowski et al. (200Cl| ). 



4. CONCLUSIONS 



We have analysed the /-mode frequencies and split- 
ting coefficients obtained from MDI medium-^ pro- 
gram for the period May 1996 to August 2000. We 
find that the mean shift in these mode frequencies 
are weakly correlated with activity indices as com- 
pared to the p-mode frequencies. In case of MWSI, 
which represents the magnetic field above 100 Gauss, 
the /-modes have the best correlation while in case of 
p-modes, the correlation is weakest. We have also es- 
timated the relative difference in the seismic radius 
as due to the change in frequencies and find that 
the change is a function of degree of the mode. If 
we consider only those modes with degrees between 
100 and 200, we find the change in solar radius from 
minimum to maximum of the solar cycle is approx- 
imately 7 km and is one order smaller than that is 



expected from the analysis of MDI data by |Dziem 
Ipowski et al. (1998|) . If we use the high degree modes 



£ u p to 286 which is ap proximately the same as used 



by Dziembowski et al. , we find a significantly higher 
value. We also report here a 1 year periodicity in 
both the /-mode frequencies and seismic radius, al- 
hough it is not clear if these are associated with the 
solar cycles or reflects the systematic errors in the 
data. 
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